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ABSTRACT 

One of the most energetic gamma-ray burst GRB 110731A was observed from optical to GeV energy range. 
Previous analysis on the prompt phase revealed similarities with the Large Area Telescope (LAT) bursts ob¬ 
served by Fermi: i) a delayed onset of the high-energy emission (>100 MeV), ii) a short-lasting bright peak 
at later times and iii) a temporally extended component from this phase and lasting hundreds of seconds. Ad¬ 
ditionally to the prompt phase, multiwavelength observations over different epochs showed that the spectral 
energy distribution was better fitted by a wind afterglow model. We present a leptonic model based on an 
early afterglow that evolves in a stellar wind of its progenitor. We apply this model to interpret the temporally 
extended LAT emission and the brightest LAT peak exhibited by the prompt phase of GRB 110731A. Addi¬ 
tionally, using the same set of parameters, we describe the multiwavelength afterglow observations. The origin 
of the temporally extended LAT, X-ray and optical flux is explained through synchrotron radiation from the for¬ 
ward shock and the brightest LAT peak is described evoking the synchrotron self-Compton emission from the 
reverse shock. The bulk Lorentz factor required in this model (T ~ 520) lies in the range of values demanded 
for most LAT-detected gamma-ray bursts. We show that the strength of the magnetic field in the reverse-shock 
region is ~ 50 times stronger than in the forward-shock region. This result suggests that for GRB 110731A, 
the central engine is likely entrained with strong magnetic fields. 

Subject headings: gamma-rays bursts: individual (GRB 110731A) — radiation mechanisms: nonthermal 


1. INTRODUCTION 

In recent years, the detection of 7 -rays and optical polar¬ 
ization in gamma-ray bursts (GRBs) has supported the idea 
that jets could be magnetized (Coburn & Boggs 2003; Boggs 
& Coburn 2003; Rutledge & Fox 2004). The jet evolution 
with magnetic content has been explored in several contexts. 
In these models, an electromagnetic component is introduced 
through the magnetization parameter (a) and defined by the 
ratio of Poynting flux (electromagnetic component) and mat¬ 
ter energy (intemal+kinetic component) (Drenkhahn 2002; 
Vlahakis & Konigl 2003a,b; Wheeler et al. 2000; Blandford 
2002; Lyutikov & Blandford 2002; Spruit & Drenkhahn 2004; 
Zhang & Yan 2011; Fraija 2014). 

Afterglow transition is one of the most interesting and least 
understood gamma-ray phases. During this phase, the rel¬ 
ativistic ejecta interacts with the surrounding matter gener¬ 
ating reverse and forward shocks. A strong short-lived re¬ 
verse shock (RS) propagates into the ejecta whereas the long- 
lasting forward shock (FS) leads to a continuous softening 
of the afterglow spectrum (Panaitescu 2007; Nakar & Pi- 
ran 2004). The dynamics of the RS in a wind and constant 
medium has been discussed by many authors (Kobayashi & 
Zhang 2003; Wu et al. 2003; Meszaros & Rees 1997a; Sari 
& Piran 1999; Fraija et al. 2012a,b; Sacahui et al. 2012; Li 
& Chevalier 2003). The RS has been invoked to explain the 
early 7 -ray, optical and/or radio flares. After the peak, no 
new electrons are injected and the material cools adiabati- 
cally, although if the central engine emits slowly moving ma¬ 
terial the RS could survive from hours to days (Genet et al. 
2007; Uhm & Beloborodov 2007). On the other hand, the 
origin of early flashes has been also discussed in the internal 
shock framework when they are nearly two orders of mag¬ 
nitude weaker than those produced by RSs (for the same to- 
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tal energy) (Meszaros & Rees 1997a, 1999; Kumar & Piran 
2000). Early observations of GRB afterglows would offer to 
clarify the question whether the early emission takes place at 
internal or external shocks. 

As has been pointed out in the literature, emission regions and 
radiative processes of high-energy (HE) photons with ener¬ 
gies >100 MeV have been fully explored. On hadronic mod¬ 
els, 7 -ray components have been explained through photo- 
hadronic interactions between HE hadrons accelerated in the 
jet and internal synchrotron photons (Asano et al. 2009; Der- 
mer et al. 2000), inelastic proton-neutron collisions (Meszaros 
& Rees 2000) and interactions of HE neutrons and photons 
out of the jet (Dermer & Atoyan 2004; Alvarez-Muniz et al. 
2004). On leptonic models, 7 -ray fluxes have been explored 
with inverse Compton (IC), synchrotron self-Compton (SSC) 
and synchrotron processes at different regions of the jet. By 
considering electrons and photons at internal and external 
shocks, IC emissions have been discussed in detail in inter¬ 
nal shocks (Papathanassiou & Meszaros 1996; Pilla & Loeb 
1998; Panaitescu & Meszaros 2000; Gupta & Zhang 2007), 
FSs (Sari et al. 1996; Totani 1998; Waxman 1997; Panaitescu 
& Meszaros 1998a; Wei & Lu 1998; Chiang & Dermer 1999; 
Dermer et al. 2000; Panaitescu & Kumar 2000) and RSs 
(Wang et al. 2001a,b; Pe’er & Wijers 2006). GeV photons 
generated from SSC emission in FS (Sari & Esin 2001; Wang 
et al. 2001a; Zhang & Meszaros 2001; Veres & Meszaros 
2012; Sacahui et al. 2012) and RS (Granot & Guetta 2003; 
Wang et al. 2001a,b; Fraija et al. 2012a) have been investi¬ 
gated separately and only in few cases synchrotron radiation 
has been examined (Liu & Wang 2011 ; He et al. 2011 ; Piran & 
Nakar 2010; Kumar & Barniol Duran 2009, 2010; Ghisellini 
et al. 2010). Additionally, a particular lepto-hadronic model 
was developed to explain the HE emissions of GRB090510 
(Razzaque 2010). 

The bright and long GRB 110731A was detected by Fermi 
and Swift observatories, and by the MOA and GROND op- 
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tical telescopes. The analysis of the prompt phase revealed 
a brightest peak in the LAT light curve starting at ~ 5.5 s 
(Ackermann & et al. 2013) and a temporally extended LAT 
component described with a power law. In addition, tempo¬ 
ral and spectral analysis in different wavelengths and epochs 
(just after the trigger time and extending for more than 800 
s) favored a wind afterglow model. Recently, assuming that 
the long-lasting LAT component could be described as syn¬ 
chrotron radiation by relativistic electrons accelerated through 
FSs and requiring that the magnetic equipartition parameter 
varies as a function of time, e# oc t~ at with 0.5 < a t < 0.4, 
Lemoine et al. (2013) interpreted that these GeV photons 
were likely produced in a region of strong e#. They ar¬ 
gued that the magnetization that permeates the blast wave of 
GRB 110731A can be described as partial decay of the micro¬ 
turbulence (Rossi & Rees 2003; Pe’er & Zhang 2006) as ob¬ 
served in particle-in-cell (PIC) simulations (Spitkovsky 2008; 
Martins et al. 2009; Haugbplle 2011; Sironi et al. 2013; Sironi 
& Spitkovsky 2011). 

In this paper, we develop a leptonic model based on early 
afterglow with variable density (stellar wind, s= 2 ) to de¬ 
scribe the temporally extended emission and the brightest 
peak present in the LAT light curve of GRB 110731A. The 
paper is arranged as follows. In Section 2 we show a lep¬ 
tonic model based on external shocks (forward and reverse) 
that evolves adiabatically in a stellar wind. In section 3 we 
apply this model to GRB 110731A as a particular case and in 
Section 4 we give a brief summary. 

2. EXTERNAL SHOCK MODEL 

As the blast wave extends out into the stellar dense wind 
of the progenitor, it starts to be decelerated leading to for¬ 
ward and reverse shocks. The evolution of the afterglow will 
mainly depend on its mass and in some cases, the emission 
processes of internal and reverse shocks which could be si¬ 
multaneously present in the light curve. In the following sub¬ 
sections we will develop the dynamics of the external shocks 
in stellar winds when Fermi-accelerated electrons are cooled 
down by synchrotron and Compton scattering emission at for¬ 
ward and reverse shocks. In addition, we will consider the 
RS in the thick- and thin-shell case. We hereafter use primes 
(unprimes) to define the quantities in a comoving (observer) 
frame and c=ft=l in natural units. The subscripts f and r re¬ 
fer throughout this paper to the forward and reverse shocks, 
respectively. 

2 . 1 . Forward Shocks 

Afterglow hydrodynamics involves a relativistic blast wave 
expanding into the medium with density 

p = Ar~ 2 with A=2L_, ( 1 ) 

47 tV w 


From eqs. (1), (2) and (3), we get the scale of deceleration 
time 


tde 


647r£ : 


-(l + z)EA- 1 T 
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Synchrotron emission — . Considering that electrons are ac¬ 
celerated to a power-law distribution N('y e )d^ e oc r y~ p d r y e 
with the electron spectral index p> 2 , and the energy den¬ 
sity (U) is equipartitioned to accelerate electrons (U e = 
e e jU = m e f j e N(j e )dj e ) and to amplify the magnetic 
field Ubj = £bjU (with Ubj = the minimum 

electron Lorentz factor and the magnetic field can be written 
as 


7e,m,/ 


(P - 2 ) mp 

(P ~ 1) m e £eJ 


r, 
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and 


Bf ^ ^^(1 + *) 1 / 2 4 /2 / r E~ 1/2 1- 1 ' 2 A , (6) 


respectively. Here, e e j (e^j) is the electron (magnetic) 
equipartition parameter and m e (m p ) is the electron (proton) 
mass. When the expanding relativistic ejecta encounter the 
stellar wind, it starts to be decelerated, then electrons are 
firstly heated and after cooled down by synchrotron emis¬ 
sion. Comparing both time scales, the deceleration time (eq. 
4) and the characteristic cooling time for synchrotron radia¬ 
tion t e ,syn ^ 3m e /(16a T ) (1 + Xf ) _1 (1 + z) p - 1 r - 3 7e -1 , 
the characteristic electron Lorentz factor can be written in the 
form 


7e,c,/ 


3m e g 4 

G T 


(1 + Xf) 1 (1 + z) 1 e B jTA l t. 


(7) 


Here gt is the Thomson cross section and the term (l+xy) is 
introduced because a once-scattered synchrotron photon gen¬ 
erally has energy larger than the electron mass in the rest 
frame of the second-scattering electrons. It is given by (Sari 
& Esin 2001) 


l+Xf 


1 + 
1 + 


V e e,f 



if 
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where 



for slow cooling, 
for fast cooling. 


( 8 ) 
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The maximum electron Lorentz factor can be calculated com¬ 
paring the acceleration t acc ~ 27r g me (1 + z) T _1 B'J 1 j e and 
cooling (t e ,syn) time scales. Hence, the maximum electron 
Lorentz factor is 


where M w is the mass loss rate and V w is the wind velocity. 
For an ultra relativistic and adiabatic blast wave, the radius 
shock (r) spreading into this density can be written in the form 

r = -24-(i _|_ 2) -1 / 2 E 1 / 2 f 1 / 2 A“i/ 2 . (2) 

Here the total energy ( E ) of the shock is constant and given by 
E = Stt/9 AT 2 r, T is the bulk Lorentz factor, £ is a constant 
parameter (Panaitescu & Meszaros 1998b), z is the redshift 
and t is the time in the observer’s frame (Dai & Lu 1998; Sari 
1997; Panaitescu & Meszaros 1998b) which is given by 

r 


'it,max,! = \7L 1/2 ( 1 + *r 1 / 4 'B 1 / 1 r- 1 / 2 E 1/1 

y lo ttgb 

xA~ 1/2 t 1/4 , ( 10 ) 

where q e is the elementary charge. From eqs. (5), (7), (10) 
and (4), the synchrotron spectral breaks computed through 
the synchrotron emission Eij = (1 + z)~ x T B'^^j for 

i=m,c and max can be written as 

E syn „ 2 16 / 5 7T 21 / 10 3^^ (p + 2) 3 / 5 (p - l ) 8 / 5 6/5 

7,a,f 5 3 / 5 r(5/6) 3 / 5 (3 p + 2) 3 / 5 (p — 2) rr?J h 

x (1 + zy 21 * e -j egj E~ 2/s 


t — (1 + 2:) 


4£ 2 r 2 ‘ 


(3) 
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The synchrotron self-absorption energy E^ y ™ f was calcu¬ 
lated through the absorption coefficient (Rybicki & Lightman 
1986) 

Q , ^ 2 22 /V 5 / 6 q e 8/3 (p + 2) (p - 1) 8/3 r 8/3 (1 ,4/3 

4 15r(5/6)(3p + 2)(p- 2) 5 / 3 m® /3 

x e “® /3 e B /3 £“ 4/3 A 8/3 t~ 4/3 E' 1 s iy 5/3 ,( 12 ) 

and the condition a e ^r/T = 1 (Wijers & Galama 1999; 
Granot et al. 1999). The maximum flux f = 

N e P lKma x /4 7 TL * 2 with peak spectral power P l/r 
CT(m e /3q e ) (1 + z)~ l T B'j can be written as 


~ ( 2 1/2 7r 3/4 m 3 g 3 (p- 1) \ / 3+2 qern^jp - 2)‘ 
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respectively. For slow-cooling regime, we get 
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(13) 

where D is the luminosity distance from the source. Follow¬ 
ing Sari et al. (1998), the observed synchrotron spectrum can 
be derived in the fast- and slow-cooling regime as 


SSC emission — . Fermi-accelerated electrons in the FS may 
scatter synchrotron photons up to higher energies E* s ? ~ 

27 \^E sy ™. From the synchrotron spectral breaks (eqs. 11) 
and break electron Lorentz factors (eqs. 5, 7 and 10), the 
spectral breaks in the Compton regime are 
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xT 5/4 t 1/4 . (22) 


respectively. Here F" 7 ^ ax f and E^ n f are given by eqs. (13) 

and (11), respectively. The transition time (t^) from fast- to 
slow-cooling spectrum is 


, syn _ gt m p (p — 2) 

0 37T 1 / 4 ra 2 (p — 1)£ 4 


(1 + z)e e ,f £bj A . 


(16) 


Using the synchrotron spectral breaks radiated (eq. 11) and 
synchrotron spectrum in the fast (eq 14) and slow (eq 15) 
cooling regime, one can obtain the light curves (LCs) as a 
function of energy (E Eyn ). We get the flux for fast-cooling 
regime 

f Asyn (psyn\-\ psyn . psyn ^ psyn 

psyn _ I 1 t^7 ) > ^ ^7 ^ ^7^4’ Q 7 ) 

\ (Err 2 . <^.p 

where and are 


Here we have calculated the self-absorption energy in the 
Compton regime through E^ f ^ & ^ s i n 

2001). From eqs. (1), (2) and (13), we get 

F 7,max,f - ( c T / m p )rp 

= + 2 ) 2 i / 2 ,4 5 ' 2 D - 2 r 1 .(23) 

In the Klein-Nishina (KN) regime, the emissivity of IC radia¬ 
tion per electron is independent of the electron energy and re¬ 
duced in comparison with the classical regime. The observed 
KN break energy is 

pKN T 

E^ f ~ —— 7 e ,c,f m e 


l + z 
, 3 mie 
gt 


(l+x f )- l (l + z)- z e- B \ f V 2 A~ l t, (24) 


where we have used eqs. (4) and (7). From the synchrotron 
spectra (eqs. 14 and 15), we get the SSC spectra for the fast¬ 
cooling regime 
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E: 
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where (£ 7 F 7 )^ >f = E^F^E^) and for the slow- 
cooling regime 
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with x/ given by eqs. ( 8 ) and (9), and {F 1 F 1 ) s ^ Ki = 
^ 7 ^f F ^ n ( £, 7 ^f) ( Sari & Esin 2001 ; Granot & Guetta 
2003). From the break photon energies (eq. 22) and syn¬ 
chrotron spectrum in the fast (eq. 25) and slow (eq. 26) cool¬ 
ing regime, one can obtain the LCs for the fast- cooling 
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and the slow-cooling regime 
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2 . 2 . Reverse Shocks 


(Usov 1992; Meszaros & Rees 1997b; Wheeler et al. 2000) 
the fireball could be endowed with ’’’primordial” magnetic 
fields. Also as the cooling Lorentz factor must be corrected, 
then 7 Z x is introduced as a correction factor for the SSC cool¬ 
ing, where x r is obtained by (Kobayashi & Zhang 2007) 


x r 


r H e e,r 



if 

if 


< 1 , 


> 1 . 


(32) 


Here r] = (7 e , c ,r/7e,m,r) 2 p is given for slow-cooling and 
r] = 1 for fast-cooling regime. 


2.2.1. Thick-shell case 

In this case, the RS becomes relativistic during its propa¬ 
gation and the ejecta is significantly decelerated. The bulk 
Lorentz factor at the shock crossing time td — Tgo is given 
by Td rsj T c . Eventually, the shock crossing time could be 
shorter than T 90 depending on the degree of magnetization of 
ejecta, defined as the ratio of Poynting flux to matter energy 
flux a = Lpf/Lkn ~ £B,r (Fan et al. 2004b; Zhang & 
Kobayashi 2005; Kobayashi & Zhang 2007). In particular, 
numerical analysis performed by Fan et al. (2004b) revealed 
that for the value of the magnetization parameter a ~ 1 , the 
shock crossing time becomes td ~ T 90 /6. 


Synchrotron emission — . Assuming that electrons are acceler¬ 
ated in the RS to a power-law distribution and the energy den¬ 
sity is equipartitioned between electrons and magnetic field, 
then the minimum electron Lorentz factor and magnetic field 
are 


For the RS, a simple analytic solution can be derived taking 
two limiting cases, thick- and thin-shell case, (Sari & Piran 
1995) by using a critical Lorentz factor (T c ) which is defined 
by 


r, = . 


8 7 T 1 / 2 £ 


(1 +z) 1/4 £ 1/4 i^ 4 T 9 “ 1/4 , (29) 


where T 90 is the duration of the GRB. For T > T c (thick shell) 
the shell is significantly decelerated by the RS 2 otherwise, 
T < T c (thin shell), the RS cannot decelerate the shell ef¬ 
fectively. Irrespective of the evolution of RS, the synchrotron 
spectral evolution between RS and FS is related by (Zhang 
et al. 2003; Kobayashi & Zhang 2007; Fan & Wei 2005; Fan 
et al. 2004a; Jin & Fan 2007; Shao & Dai 2005) 


where 
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and 72 m — 


T ’ 
(31) 


where Td is the bulk Lorentz factor at the shock crossing 
time. The previous relations tell us that including the re¬ 
scaling there is a unified description between forward and re¬ 
verse shocks, and the distinction between forward and reverse 
magnetic fields considers that in some central engine models 


2 Although bulk Lorenz factors (FS and RS) can be different at the shocked 
region, we have considered them to be similar. 
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From the characteristic cooling time of synchrotron radiation 
and dynamical time scale, the characteristic electron Lorentz 
factor can be written as 
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By considering 7 e?a ~ 7 e?m (Sari & Esin 2001) and from 
eq. (31), we re-scale the synchrotron self-absorption energy 
between FS and RS as Eppp ~ 7 Z 2 e K~ 1/5 Kpp E° v £ f . Addi- 
tionally, from eqs. (11), (31) and (30), we get the synchrotron 
spectral breaks 
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and the peak flux 


F Syn ~ y/ ^ me<JT (i I z )2 1/2 p -1 jl /2 jj-2 
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Synchrotron LCs are derived by Kobayashi (2000). Flux 
increases proportionally to ^ t 1 / 2 reaching a peak at 

FZ n e ak,r ~ (^777) “ 1/2 and after it starts de- 

creasing as ~ t~ 3 (Kobayashi & Zhang 2003). 


SSC emission — . Accelerated electrons can upscatter photons 
from low to high energies as 
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where k = 4(p - l)/(p - 2 ) and r = CT 7 4 ^J e) is the optical 
depth of the shell. Here 7V e is the number of radiating elec¬ 
trons . From eqs. (38), (44), (33) and (47), we get the break 
SSC energies 


JjjSSC ^ 
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xD~ 2 ET 9 ~ 0 2 ,( 39) 


and the break energy at the KN regime is 

(i+^r 1 (l+Xr+xl)- 1 e- B ] r ET~ 2 A - 2 . 

(40) 

At the shock crossing time (td), the SSC flux reaches the peak 

7,peak,! ~ (E%/E% c c ,r)- 1/2 i^,r(Kobayashi & Zhang 

2003) at: 


rpssc 

-*7,peak,r 




X € e 


-7/2 p -6 a- 
’ e S,r 1 ^ 


6 d~ 2 et~ 0 1/2 {e™?} 


- 1/2 


(41) 


2.2.2. Thin-shell case 


In the thin-shell case, the RS cannot decelerate the shell 
effectively. The deceleration time and the minimum electron 
Lorentz factor are in the form 


and 



7e,m,r — ^ 


(i + ^)^^- 1 r- 4 , 

(42) 

fp-2\m p 
e ’ r \p-l)m e ’ 

(43) 


respectively. Here the bulk Lorentz factor at the shock is T 

r d < r c . 


Synchrotron emission — . Performing a similar analysis of the 
thick-shell case, from eqs. (11), (31), (30) and (43), we get 
that the synchrotron break energies and maximum flux can be 
written as 
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(44) 


1 S 1 SC emission — . In a like manner to the thick-shell case, from 
eqs. (38), (44), (33) and (47), we derive the break energies of 
SSC emission 
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SSC LCs can be analytically derived from Chevalier & Li 
(2000). For t < td, we take into account that: i) the num¬ 
ber of radiating electrons N e and the spherical radius in the 
shocked shell region increase with time as N e ~t and r ~ t, 
ii) the maximum flux of synchrotron is independent of time 
-^ 7 ,max 00 an d iii) the SSC cooling break energy E* s ° ~ 
7 cEfv™ increases as ^ t 3 , then the SSC flux increases as 
F ssc ~ F ssc1 / 2 F ssc ~ F ssc1 / 2 —^F s yn ^ y-i/2 p or / \ 

-4- ij r x 1 / m fiy 0 / r ^.2 0 / ma y u • 1 u ^ 


7,c 


7,c 


7 , max 


t d , the SSC flux decreases as F° sc - ^“-(p- 1 )/ 2 ]?^ ~ 

inssc - ^ 5 - iVe rpsyn j . -£±1 where K S5C ~ 'v 2 E syn de- 

- / - / 7 ,m r 2 - r 7 ,max 6 ’ Wlieie 7m^7,m 

creases as ~ f _1 . It is valuable to note that the decay index of 


the emission for t > td might be higher than ^ > due to the 

angular time delay effect (Kobayashi & Zhang 2003). 


The break energy at the KN regime is given by eq. (40). 

3. APPLICATION: GRB 110731A 

GRB 110731A was detected on 2011 July 31 by both 
instruments aboard Fermi; Gamma-Ray Burst Monitor 
(GBM) and Large Area Telescope (LAT) (Ackermann & 
et al. 2013), the three instruments aboard Swift; BAT, 
XRT and UVOT (Oates et al. 2011) and ground-based 
observatories; Microlensing Observations in Astrophysics 
(MAO) telescope (Tanvir et al. 2011) and Gamma-ray Burst 
Optical/Near-Infrared Detector (GROND). 

LAT localized GRB 110731A with coordinates R. 
A _i8^41 m 00 s and dec.=-28 o 31’00” (J2000), with a 
68 % confidence error radius of 0.2°. Swift/BAT perceived 
immediately this burst after the detection by both instruments 
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of Fermi, whereas RXT and UVOT began observations 56s 
after the BAT trigger. UVOT swiftly determined the afterglow 
position as R.A. = 18' l 42 m 00 s .99 and dec.=-28°32’13”.8 
(J2000), with a 90% confidence. The lack of observation 
in the UV filters is consistent with the measured redshift 
z=2.83 (Tanvir et al. 2011). MAO observations started 3.3 
minutes after the Swift trigger for GRB 110731 A. Using a 61 
cm Boiler & Chivens telescope at the Mt. John University 
Observatory in New Zealand, I- and V-band images were 
collected 105 minutes after the trigger (Tristram et al. 2011). 
Finally, GROND mounted on the 2.2 m MPG/ESO telescope 
at La Silla Observatory, Chile, imaged GRB 110731A for 
2.74 days after the trigger (Greiner et al. 2008). 

By considering the typical values of the stellar wind 
(A = A* x (5.0 x 10 11 ) gram/cm with A * = 0.1; Chevalier 
& Li (2000), the parameter £ = 0.56; Panaitescu & Meszaros 
(1998b)) and those inferred by observations: redshift 
z = 2.83 (Tanvir et al. 2011), total energy E ~ 10 54 erg and 
duration of GRB T 90 = 7.3 s, we will apply the leptonic model 
developed in this work to interpret the LAT LC observations. 
Taking into account the fact that the peak of the flux density 
was present at the end of the prompt phase; in the interval 
[5.47 s, 5.67 s], and that after the LAT flux decays smoothly 
during the whole temporally extended emission, we constrain 
the Lorentz factor T ~ 520 > T c so that the deceleration 
time takes place at 

tdec - 5 . 55 s^Lti:^ Er-,4 Tp r2 , (46) 

and the RS evolves in the thick-shell case with a critical 
Lorentz factor 

472 . 5 (UA) ' ' - (47 > 

We plot the synchrotron and SSC spectral breaks of the FS 
and RS as a function of equipartition parameters (e B j^ 
and e e j( r )), considering the typical values of the magnetic 
(10 -5 < CBj(r) < 1) an d electron ( e e j( r )= 0.5, 0.1 and 
0.05) parameters (Santana et al. 2014), as shown in fig 1. We 
describe this figure as follows. 

Synchrotron spectral breaks from FS (right-hand panel). — From 
this panel, we can see that the characteristic energy is an in¬ 
creasing function of e B j and e e j, the cooling energy is a de¬ 
creasing function of e B j and the self-absorption energy is an 
increasing function of e B j and a decreasing function of e e j. 
The characteristic energy lies in the ranges: 5.3 x 10 2 eV < 
E^Zj < 2.8 x 10 5 eV, 3.5 x 10 3 eV < E sy ^ f < 6.8 x 10 5 eV 
and 6.1 x 10 4 eV < E sy ^ f < 1.9 x 10 7 eV for e eJ = 0.05, 
0.1 and 0.5, respectively. Additionally, the cooling energy 
lies in the range 8.9 x 10 -6 eV < E* v £ f < 2.1 x 10 eV 
and the synchrotron self-absorption energy in the 
ranges 5.2 x 10“ 3 eV < E* v ” jf < 4.6 x 10 -2 eV and 

5.1 x 10 -4 eV < E° y 2 f < 6.3 x 10“ 3 eV for e eJ = 0.05 
and 0.5, respectively. Furthermore, one can see that self¬ 
absorption energy is in the weak self-absorption regime 
( E * y aj < ESy cj) for £BJ < 0.003 (0.007) and e eJ = 0.05 
(0.5), otherwise the synchrotron spectrum would be in the 
strong absorption regime. 


RS have the same behavior as FS although at lower energy 
ranges. The characteristic synchrotron energy is in the ranges: 
0.01 eV < E°% tr < 5.1 eV, 0.08 eV < E* y ^ r < 1.1 x 10 2 eV 
and 1.2eV <’ E°% tr < 6.2 x 10 2 eV for e’, r =0.05, 0.1 and 
0.5, respectively. The cooling energy is in the range 1.8 x 
10 _6 eV < E sy Jf^ r < 46.2 eV and the self-absorption energy 
in the ranges 3.8 x 10 _7 eV < E syr ^ r < 8.1 x 10 _7 eV and 
7.3 x 10 “ 9 eV < E*y™ r < 6.8 x lO” 8 eV for e e , r = 0.05 and 
0.5, respectively. It is important to say that synchrotron self¬ 
absorption is in the weak absorption regime (E^ y ™ r < E^ y ™ r ) 
for any value of equipartition parameters considered here. 
However, as E*v™ r oc A 11 / 5 and E sy ^ r oc A~ 2 any signifi¬ 
cant increase of wind density, the absorption energy would be¬ 
come higher than cooling energy (E* y ™^ > E^ y ^), then the 
synchrotron spectrum could change to the strong absorption 
regime. In this case, heating of low-energy electrons due to 
synchrotron absorption leads to pile-up of electrons, and there 
appears a thermal component besides the power-law spectrum 
(Gao etal. 2013). 

SSC spectral breaks from RS (panel below). — SSC spectral 
breaks have similar behavior as in the previous cases. The 
characteristic SSC energy lies in the range 52.3 eV < 
Ey r m r < 1.12 x 10 4 eV, 7.8 x 10 3 eV < E^ n r < 2.3 x 
10 5 eV and 6.2 x 10 5 eV < Ey^ r < 1.17 x’l0 8 eV for 
e e r =0.05, 0.1 and 0.5, respectively while the cooling energy 
lies in the range 8.2 x 10 -1 ° eV < E^ s f r < 2.1 x 10 10 eV. 
To obtain the values of parameters e B j( r ) an< i e ej(r ) that 
reproduce the multi wavelength LC observations, we use the 
method of Chi-square (% 2 ) minimization (Brun & Rademak- 
ers 1997). We describe LAT flux observations by synchrotron 
radiation from FS and SSC emission from RS; the whole tem¬ 
porally extended emission using synchrotron LCs in the fast¬ 
cooling regime (eq. 17) and the brightest peak at 5.5 s with 
SSC LCs in the the fast-cooling regime (eq. 41), for high- 
energy electrons radiating at ~ 100 MeV. Following the anal¬ 
ysis showed in Ackermann & et al. (2013), we fit the X-ray 
(t < 4.6 ks) flux with the synchrotron LC in the fast-cooling 
regime (eq. 17) at t= 100 s for electrons radiating at E sy ^ = 
2 keV, and optical and X-ray (t > 4.6 ks) fluxes with the LC 
in slow-cooling regime (eq. 20) at t=600 s and 4000 s for 
E SyT f ~ 10 an d 0-7 keV, respectively. 

We plot the values of parameters e B j( r ) and e e j( r ) for 
p=2.15, 2.2, 2.25 and 2.3 (see figs. 2 and 3) that reproduce 
the multiwavelength LC observations (see fig. 4). As shown 
in fig. 2 the areas in yellow (green) colors show the set of 
parameters that describes the temporally extended LAT (op¬ 
tical) flux and in blue (red) ones show those parameters that 
describe the X-ray for t < 4.6 ks (t > 4.6 ks). The areas in 
brown colors show the set of parameters that reproduce more 
than one flux at the same time. For instance, as shown in 
the right-hand panel below (p=2.25), the equipartition param¬ 
eters in the narrow strip between 0.38 < e e j < 0.52 and 
IQ - 4 5 < eB j < 10 -4 would reproduce the temporally ex¬ 
tended LAT, X-ray and optical fluxes. From fig. 3, one can 
see the set of parameters that describe the brightest LAT peak. 
Also it can be seen that as the electron spectral index increases 
the set of parameters is shifted to the right. 

Given the values e B j = 10 -415 and e e j = 0.40 for p=2.25, 
from eqs. ( 11 ) and (16) we get that the synchrotron spectral 
breaks from FS are 


Synchrotron spectral breaks from RS (left-hand panel ).— From 
this panel, we can see that synchrotron spectral breaks from 


^W-5.56 x 10 -4 eV 


1 + z 
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77,36.94 GeV t 


-3/4 


x£ s ‘f«; /s 


1 + ^ \ 7^1/4 a — 1/4 .—1/4 

^54 ^*,-1 r i 


(48) 


and the transition time from fast- to slow-cooling regime is 

j_syn _ 223.465 ^^ e e J,-0A CBJ,-4.15 ^*,-1 • (49) 

It is important to highlight that the maximum photon energy 
achieved by synchrotron radiation is f — 36.94 GeV 

for t=10s and L^ ax?f ~ 20.77 GeV for t=100s. 

From eq. (26) we get that the SSC scattering break energies 
are 


11.66 TeV 


f 1 + Z \ ^4 ^ 1/2 




e e,/,-0.4 e S,/,-4.15 

X A^J_i F /54 f 2 2 
-3 / \ -4’ 


^162.8keV^l7j (^j e-:f _, 15 

xA:^E 54 tl 

F°% ax.f ^ 0.41 Jy 4 / ,/,-4.i5 V-1 7^54 • (50) 

The break energy at the KN regime is = 42.33 GeV. 
From the SSC LCs (eqs. 27 and 28), one can see that al¬ 
though the temporal power index of LC for f < E^ sc < 
f would be in agreement with the temporally extended 
LAT flux, the energy range would not. The detection of the 
photon with energy 3.4 GeV at ~ 436 s could be explained by 
synchrotron radiation as well as Compton scattering emission. 
Considering e e?r = e e j = 0.4, we obtain €bj = 0.28 for 
p=2.25. From eqs. (44) and (39) we get that the synchrotron 
and SSC break energies are 


\ “ 7/5 

7a"r ^ 4.28 X IQ" 8 e V ( ~E~ I ^ejr,—0.4 44,-0.55 ^2.72 


e;I /5 (44' 


\ 2 / 5 


7.3s ) 


\ _1/2 

E%Z,r - 128.94 eV ( ^ ) e 2 ,,,_ 0 . 4 44 - 0.55 41.72 


(is) 


- 1/2 


E* v £ t ~ 0.93 x 10 eV x 


\ _3/2 / 2 
—5 / 1 T z \ / 1 + 37 T 37 


x f~ 3/2 A~ 2 E 1 ' 2 ( — 

XC B,r,- 0.55 ^*,-1 -^54 (735 


1/2 


-F 75max ,r 4.31 x 10 Jy ( ^ IV 72 ^J,-i 


Tqq \ 


xD^E 54 ( Eg) , (51) 


and 


77 ~ 103.55 MeV ( ) 7,-o.4 44.-o.55 U 72 

y A 3 / 2 E~ X 
X/1 *,-l ^54 5 

. -3/2 / 2 \ - 4 

T? SSC r Qa ^ 1 n -3 -t r / I H - ^ \ / 1 + X r + X r 

L/ 7;C;I . ~ 5.86 x 10 eV [ —— 


-7/2 p — 6 a —6 p3/2 \ 

3,r,-0.55 1 2.72 ^*,-1 ^54 l ^ 3 ^ ) ’ 


xe 


7W ~ 1.42 x 10 2 Jy ( 1^ I 44-0.55 r 2 : 2 72 

-2 (t» o_ N _2 


xD 2 7B h ^ ] , (52) 


respectively. The break energy at the KN regime is E™ = 
52.71 GeV. 

4. CONCLUSIONS 

We have presented a leptonic model based on the evolution 
of an early afterglow in the stellar wind. We apply this 
model to describe the temporally extended emission and 
the brightest peak present in the LAT light curve of GRB 
110731A, although additionally we fit the X-ray and optical 
light curve afterglows. 

We consider that the ejecta propagating in the stellar wind 
is early decelerated at ^ 5.5 s and the RS evolves in the 
thick shell regime. Taking into account the values of redshift 
2 = 2.83 (Tanvir et al. 2011), energy E ~ 10 54 erg, duration 
of GRB Xgo= 7.3 s (Ackermann & et al. 2013; Oates et al. 
2011; Tanvir et al. 2011; Tristram et al. 2011) and the stellar 
wind A = 5.0 x 10 10 gram/cm (Chevalier & Li 2000), we 
get that the value of the Lorentz factor is T ~ 520. 

We plot the SSC and synchrotron spectral breaks from FS 
and RS as a function of equipartition parameters, as shown 
in fig 1. We found that the synchrotron self-absorption 
energies from FS and RS are in the weak self-absorption 
regime for cbj < 0.003 (0.007) and e e j = 0.05 (0.5) and for 
10 -5 < €b, v < 1 and 0.05 < € e r < 0.5 (considered here), 
respectively. It is important to say that for 7 e , a > 7 e ,m (not 
considered in this work) the synchrotron spectrum would be 
in the strong absorption regime (E*v ™ f > ). In this case a 

thermal peak due to pile-up of electrons would appear around 
E syn ~ 6 23eVe 2//7 e 1//14 A 2 ^ 7 E 1 ^ 14 (2m ) -11 / 14 

^ 7 ,a,f O.Z4)e\/ fc e ,r,-0.4 6 B,r,-0.55 ^*,-1 ^54 V 7.3s / 

in the synchrotron spectrum of the RS, modifying the broken 
power-law spectrum (Kobayashi et al. 2004; Gao et al. 2013). 
To find the equipartition parameters e B j(r) an d £ e ,f(r) we fit 
the multi wavelength afterglow LCs (see fig. 4); the brightest 
LAT peak by SSC emission from RS and the extended tem¬ 
porally emissions (LAT, X-ray and optical) by synchrotron 
radiation from the FS. We plot the set of values ^Bj(r) an d 
e e j( r ) that describes these observations, as shown in figs. 2 
and 3. From the FS, we chose the values cbj = 10 -415 and 
e e j = 0.4 for p=2.25, then we get the synchrotron spectral 
breaks: E^ f ~ 5.56 x 10- 4 eV, E s ^ ~ 77.45 keV, 
£^ f ~ 0.30 eV and £^ ax?f ~ 36.94 GeV, and SSC spec¬ 
tral breaks: 7^m,f — 11.66 TeV and — 162.8 keV. 

From the RS, we consider the values e e?r = e e j = 0.4 
and €b, v = 0.28, then we obtain the synchrotron spectral 
breaks:’ E^ r ~ 4.28 x l(r 8 eV, E^ r ~ 128.94eV 
and 77 — 0-93 x 10 _ 5 eV, and SSC spectral breaks: 
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>r ~ 103.55 MeV and E^ r ~ 5.86 x 10“ 3 eV. 

The current model accounts for the main temporal and 
spectral characteristics of GRB 110731A, having six free 
parameters (bulk Lorentz factor, density of the stellar wind, 
electron and magnetic equipartition parameters). From the 
parameters required and observed (T, A, E , X 90 , p and z), 
and the values obtained after describing the afterglow LC 
( e e j = e e , r and cbj 7 ^ es, r ), we can see that the difference 
between the SSC and synchrotron spectral breaks achieved 
in forward and RSs can be explained through the energy 
distribution given to magnetic field. Comparing the magnetic 
equipartition parameters that best describe the emission at 
forward and RSs, we can see that magnetic fields in both 
shocks are related by Bf ~ 2 x 10 -2 L> r . The previous 
result as found in GRB 990123, GRB 021211, GRB 980923 
and GRB 090926A (Zhang et al. 2003; Fraija et al. 2012a; 
Sacahui et al. 2012) illustrates that the magnetic field in the 
reverse-shock region is stronger (^50 times) than in the 
forward-shock region which indicates that the ejecta is mag¬ 
netized. Unlike the FS emission that continue later at lower 
energies, the RS is shown as a single peak. After the RS has 
crossed the ejecta, there are no freshly accelerated electrons 
injected and the emission drops sharply. The magnetization 
of the blast wave modifies the temporal characteristics of the 
brightest LAT peak; it becomes much shorter in duration, for 
GRB 110731A the peak duration is ~ Is (Fan et al. 2004b). 
Some authors have claimed that the GeV emission detected 
by LAT during the prompt phase (before Tg 0 ) has an internal 
origin similar to its MeV counterpart (Maxham et al. 2011; 
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However, in particular for GRB 110731A, we describe 
the multiwavelength data with synchrotron radiation as 
originating from FS in stellar medium starting at ^ 5.5 
s. Then, it is natural to think that the brightest LAT peak 
around the afterglow onset time comes from the RS as was 
explained in this work. We note that bursts with sub-TeV 
photons at hundreds of seconds from the prompt phase could 
be described by SSC emission from FS and be candidates to 
be detected by TeV 7 -ray observatories as the High Altitude 
Water Cherenkov observatory (HAWC) (Abeysekara & et al. 
2012, 2014). 

In summary, we model not only the LAT light curve (the 
long-lasting GeV emission extended up to 853 s by syn¬ 
chrotron emission from FS and the brightest peak by SSC 
emission from RS) but also explain the multiwavelength 
afterglow observations in GRB 110731A using the leptonic 
model based on an early afterglow which evolves in a stellar 
wind. The bulk Lorentz factor required in this model is 
T ~ 520 and the ejecta must be magnetized. In this model 
the onset of the HE emission is delayed because it is emitted 
from external shocks. 
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Fig. 1.— Break photon energies of synchrotron radiation from the forward (left-hand figure above) and reverse (right-hand figure above) shocks, and SSC 
emission from the RS (figure below) as a function of magnetic equipartition parameter (e B j/ r ) for e e j / r = 0.5, 0.1, 0.05 and 0.01. 
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Fig. 2.— Values of equipartition parameters (cbj and e e j) that reproduce the temporally extended LAT, X-ray and optical emission through synchrotron 
radiation from FS. 
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Fig. 3.— Values of equipartition parameters (e^,r and e e ,r) that describe the brightest LAT peak through SSC emission from RS. 
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Fig. 4.— Fits of the multiwavelength LCs of GRB 110731A observation with our model. We use the RS in the thick-shell regime to describe the brightest peak 
of the LAT flux (continuos line) and the FS to explain the temporally extended LAT, X-ray and optical emissions (dashed lines). 















